Branching rate expansion around annihilating random walks 
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We present some exact results for Branching and Annihilating Random Walks. We compute the 
non-universal threshold value of the annihilation rate for having a phase transition in the simplest 
reaction-diffusion system belonging to the Directed Percolation universality class. Also, we show that 
the accepted scenario for the appearance of a phase transition in the Parity Conserving universality 
class must be improved. In order to obtain these results we perform an expansion in the branching 
rate around Pure Annihilation, a theory without branching. This expansion is possible because we 
manage to solve Pure Annihilation exactly in any dimension. 

PACS numbers: 05.10.Cc 64.60.De 64.60.ae 82.20.-w 



Active-to-absorbing phase transitions represent one of 
the simplest cases where genuine non equilibrium behav- 
ior is attained. In this context, much work has been de- 
voted to the study of Branching and Annihilating Ran- 
dom Walks (BARW) [l|-|5|, systems composed of parti- 
cles of a single species A, that diffuse in a d-dimensional 
space, and that can suffer both annihilation and branch- 
ing (i.e. offspring creation) processes. BARW are not 
only of direct physical interest, but also present a con- 
ceptually simple class of out of equilibrium systems. 

Due to universality, it is generally enough to consider 
the simplest possible reactions, which allow for BARW 
to be classified into sub-classes il, 2]. Pure Annihilation 
(PA), the theory without branching where the only re- 
action is 2A — > 0, constitutes a good starting point in 
order to study properties of BARW. In the long time 
limit the PA system approaches the empty state, where 
all density correlation functions are zero. The response 
functions, however, are non trivial, and are governed in 
the Infrared (IR) (that is, for momenta and frequencies 
smaller than the scale set by A) by a non-Gaussian RG 
fixed point. If we add the simplest branching reaction 
A — > 2A, the resulting BARW system is known to be 
in the Directed Percolation universality class [Gj , and we 
call it BARW-DP. If instead we choose to add the reac- 
tion A — > 3j4, which preserves the parity of the number 
of particles, the system belongs to the Parity Conserv- 
ing (PC) universality class (more appropriately known as 
Generalized Voter [7| ) whenever a phase transition takes 
place. We call this system BARW-PC [§]. 

Both universality classes have been studied using per- 
turbative RG methods [2|, by expanding the model 
around the Gaussian reaction-less RG fixed point. In 
the case of DP, the authors of [31 do not find a phase 
transition in dimensions d > 2. As for PC, at 1- and 2- 
loop orders an (upper) critical dimension dc > 1 is found, 
with the PA fixed point being unstable for d > dc, so 
that the branching perturbation is always relevant, there 
is no absorbing phase, and hence no phase transition. 



Conversely, for d < dc the PA fixed point is fully attrac- 
tive and an absorbing phase exists at small branching a, 
whereas at larger a the system is in its active phase. A 
phase transition must therefore occur at finite a, which 
is confirmed numerically in d = 1. In [10|, a 1-loop anal- 
ysis of BARW with Levy-flight dynamics was shown to 
be compatible with this scenario. Finally, Non Pertur- 
bative Renormalization Group (NPRG) studies are also 



consistent with it [11 1 



In this letter, we revisit these conclusions, proving 
that: (i) in DP there is a phase transition for all d, in 
agreement with what is found in Montecarlo and NPRG 
studies 12|, |l3|, and (ii) in PC, the PA fixed point re- 
mains unstable in the cr-direction for all relevant dimen- 
sions d > 1, and therefore there are aspects of the PC 
transition that are still to be fully understood. 

Our results rely on an expansion in a around the PA 
theory (contrary to the usual perturbative expansion, 
which is performed around the Gaussian theory). This 
expansion is highly non trivial since it requires to solve 
exactly PA, which we do by deriving closed equations for 
all its response functions. As far as we know, such an 
expansion around a non-Gaussian model has never been 
performed before for out of equilibrium systems. Since 
our approach is valid for any value of the annihilation 
rate A, we obtain exact results at small a. 

Being based on an exact solution of PA, our approach 
also allows us to exactly compute non-universal quanti- 
ties. We choose to calculate, for the BARW-DP model, 
the non-universal threshold value Xth of the annihilation 
rate, at which a phase transition occurs at vanishing a. 

Pure Annihilation. We briefly recall some technical 
aspects. In order to study reaction-diffusion processes, a 
fleld theory can be constructed in a standard way by us- 
ing the Doi-Peliti formalism [ij]. As a result, one obtains 
the generating functional of the correlation and response 
functions 



Z[J, J] ^ / V(l3V(l)exp -S[(j), (j)]+ / J0 + J(/. (1) 



with an appropriate action S[4>,(j}], that captures the mi- 
croscopic reactions. Here we have introduced the nota- 
tion, to be used throughout, x — (x, t) (and p ~ (p, v) in 

Fourier space), J^ = J d'^xdt andJ^^J^^. 
In the case of PA the only reaction is 2 A — > and 



^^■^[0,0] 



arrives at the known equation IT], [l 



cj){dt - D\/')<j) + Xi(b' - IW ] (2) 



where D is the diffusion constant and where we have 
ignored initial conditions, because we are only interested 
in the steady state. This theory has a strong resemblance 
with equilibrium 0^ theory, but is further constrained by 
causality properties [l| . 

All the information of a system is encoded in the ver- 
tices r'"'™' of the theory (also known as the 1-Particle 
Irreducible (IPI) functions, with n incoming legs and m 
outgoing response legs), related to the connected corre- 
lation and response functions by a Legendre transform 
[15| . We now present an identity allowing us to obtain a 
closed equation for any p'"''") in PA. It can be most con- 
veniently written at the diagrammatic level: any Feyn- 
man diagram contributing to P^"'™) that includes at least 
one loop must begin by a 4-legs bare vertex. 




n — 2 legs 



FIG. 1. Generic form of a diagram for r'"'™' including at 
least one loop. Left hand side: diagrammatic representation 
of a generic r'"''"' vertex. Right hand side: general structure 
for r'"'™' in PA, the black blob is a connected and amputated 
Green function that has to comply with some requisites, see 
text. 

In Fig. [T] we show the general structure of these dia- 
grams. The black blob denotes a sub-diagram that is con- 
strained by the condition that the full diagram must be 
IPI. In particular, it means that this sub-diagram must 
be connected (and with amputated external legs). Now, 
any connected diagram contributing to this blob has a 
unique tree decomposition in terms of IPI sub-diagrams 
having at most n incoming and m outgoing legs [15| . By 
summing all possible diagrams and permutations com- 
patible with the IPI structure of the full diagram, we ob- 
tain a closed equation that relates any p("'™) with other 
pC'S) ^[f^iy g^ lower number of legs. A non-perturbative 
proof (not based on an all-order analysis) of this general 
property will be given elsewhere J16l | . 

In order to be concrete, let us analyze the identity 
given in Fig. [l]for the simplest vertices. For P^^'^) this 
yields a well-known property: there is no correction to 
p(i.i) in PA, since there is no diagram such as the one in 
Fig. [T]with a single incoming leg. Concerning P*^^^^', one 



r('^'nPi,P2,Pi,P2) = 4A-2A fciq) 

Jq 

X G{pi +P2 -g)p(^'^^((7,pi +P2 -q,Pi,P2) (3) 

with G{q) — [P'^^^'(— g)]^^ the propagator of the theory. 
The solution of ([3]) is of the form 



r'^'^'^\pi,P2,pi,P2) = 4:i{pi +P2) 



with 



l{p) = A 1 + 2A / G{q)Gip - q) 



(4) 



(5) 



This result can simply be seen as stemming from a ge- 
ometric sum over bubbles. A similar relation is found 
for P*^^'^^. The relation presented in Fig. [T]is a general- 
ization of these known results to any vertex function. It 
enables us to study BARW by means of a perturbative 
expansion in the branching rate a around PA. 

BARW-PC. Let us begin by considering the Parity 



Conserving universality class: 2A — > and A 
The corresponding microscopic action reads 



3A. 



S 



PC\ 



^{dt~DV'^)(j)+\{4>'^~l)(j)'^+a{l~^'' 



This action is symmetric under (p — > — </), cj) -^► —0, which 
implies the constraint P^"'"*) = if (n -f m) is odd. 

The PC model is known to present an active-to- 
absorbing phase transition in rf = 1, generally believed 
to be related to a change of stability of the PA fixed 
oint in a dimension dc with 2 > dc > 1- Perturbatively 
and also within the Local Potential Approximation 
(LPA) of the NPRG III , this change of stability occurs 
in the following way: On the one hand, for dimensions 
close to two, canonical power counting arguments show 
that the PA fixed point is unstable in the cr-direction, 
which implies that the system is in its active phase for 
all (7 > 0. On the other hand, at 1- and 2-loop orders 
an (upper) critical dimension dc > 1 is found such that 
for d < dc the coupling a becomes irrelevant around the 
PA fixed point, which therefore becomes fully attractive, 
thus showing the existence of an absorbing phase at small 
a. This change of stability occurs because a new fixed 
point F^^ crosses the PA fixed point at dc and in this 
dimension they both change their stability. Below dc, 
this new fixed point is in the physically relevant quad- 
rant A > 0, (T > 0, has one unstable direction, and is 
thus associated with the PC phase transition. We now 
reconsider this scenario. 

We can perform a systematic expansion in a while 
keeping A finite. Within this formalism we can reana- 
lyze the stability of the IR PA fixed point in the presence 



of the PC creation reaction, A — > 3A, that we can deter- 
mine exactly since our analysis is exact at small a. The 
relevance of this coupling can be obtained from the flow 
of r^^''^), since this function is of order a. 
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FIG. 2. Closed equation for r'^-^^ in BARW-PC. 

At first order in a, we can write the equation for P^^'^) 
in the form shown diagrammatically in Fig [21 whose 
structure implies the following functional form for pf^--^) 

r^^'^'>{p,Pi,P2,P3) = -2a{p,pi) - 2a{p,p2) - 2a{p,p3). 

(6) 
The quantity that interest us is do- , the scaling dimension 
of a in the IR limit, which can be extracted from the 
behavior of d-{p) = a{0, —p)/l{p) 

^(p)^|p|'~'", Kp)-\v?-'' for ^,|pp«A^. 
The equation for a obtained from Fig[2]and Eq. ([6]) reads 



a{p) 



a 
A 



A I G{q)G{p~qMq)l{q). (7) 



Using this exact expression and expanding in e = 2 — d 
we easily recover the 1- and 2-loop results [2, UM, ■ 

In order to get a result for da it is convenient to get 
rid of the bare reaction rates, as we are only interested 
in the universal IR scaling behavior. We subtract to ([7]) 
its value at zero a{p — 0), which is zero for d < 2, given 
that we expect d^ < d. Our results will later confirm 
this. We must also take into account the expected scaling 
invariance. We define the scaling function cr(£') 



'^(p,'^)-ipr"''"'^('> 



^- (8) 



Observe that we are performing a perturbation around 
the PA fixed point, whose anomalous dimensions are 
zero, and the natural scaling variable is accordingly 

The ensuing equation for a is still too complicated to 
be solved analytically, and requires a numerical solution. 
A convenient way to do this is to make an expansion in 
u — cos (p, q), which turns out to be rapidly convergent. 
We then proceed as follows: at each order in u, we ad- 
just da- at a given value of d, by numerically iterating this 
equation in order to reach a fixed functional form for a{v) 
in a lattice of N^, points with a resolution 5v. We have 
checked the convergence in u and in the numerical param- 
eters 5i^ and iV^, used for the computation of integrals. 
This procedure gives always a converged scaling func- 
tion a'{i>), which confirms a posteriori the scaling form 



ansatz Q. We observe a{v') to be a non-trivial function 
of its argument |16j . which may explain the qualitative 
difference between our results and previous approximate 
results. Observe that LPA and 1-loop analysis are based 
on a P^^'-^) vertex without dependence on frequency and 
momentum. 



d 0.5<-*- " ■" 




FIG. 3. Results for dcr, showing there is no change in the RG 
relevance for the branching rate a. 



This procedure allows us to find the value of da as a 
function of d, the results of which are plotted, together 
with previous perturbative results, in Fig. [31 There one 
sees that, even if da gets smaller when d decreases, it 
remains always positive, which is in contradiction with 
the usual picture for the PC transition. Notice that this 
result does not rule out the possibility of a transition in 
d — 1. We propose the following scenario that reconciles 
all existing results. In a dimension between 1 and 2 two 
fixed points appear at positive A and a, the one with 
the smaller value of a being fully attractive and govern- 
ing the absorbing phase, while the other is once unstable 
and is thus associated with the PC transition. In this sce- 
nario, the absorbing phase must have a behavior different 
from PA. This can be studied either by using Montecarlo 
methods or higher order NPRG equations. 

The previous result is surprising because another exact 
analysis [19[ in d = 1, which seems to indicate that a 
is irrelevant with respect to the IR PA fixed point, in 
contradiction with our conclusions. We can explain this 
difference observing that the model used in [19| is defined 
with A = oo (and indeed presents no phase transition at 
all for whatever value of a). Now, the IR limit of the 
theory corresponds to i^, |pp <^ X^/i^-d) ^ but this does 
not allow us to take ct = when compared to A. Looking 
at Eq. (O, A = oo implies (t = 0, so that the relevant 
direction corresponding to a is no longer accessible by 
studying ct as a perturbation. This is true for all d. In 
particular, for d ^ 2, the form of the relevant direction 
can be calculated perturbativelly and the result is at 1- 



loop 



d 



6 



a{p) ^ ^l\p) 



(9) 



One observes then that in the hniit A — )■ oo the relevant 
direction is eliminated artificially even at d ~ 2. Thus, 
the results of 19] do not apply to BARW-PC at finite 
values of the reaction rate A, the system in which we 
are interested in. It also shows that when A is large, a 
cross-over must occur and for a long transcient the PA 
behaviour will show up. Montecarlo studies of the low 
branching regime of this system have until now, as far 
as we know, also been mostly made in the limit A — > oo 
20|,|21|. They are compatible with the standard scenario, 
but within the criticims previously pointed out. 

BARW-DP. Let us now consider the simplest 

BARW-DP model: 2A ^ and A A 2A, whose mi- 
croscopic action reads (2| 



gDP^ 



(4>{dt-DV'')(j) + \4,{4> + 2)(f' 



I) 



Notice that in this equation we have performed, as is usu- 
ally done, a shift in the response fields, (j){x) = 1-1- (^(a:). 
The case <t = corresponds to PA, now written in terms 
of the shifted (j) field. This version of PA can be solved 
following the same ideas as previously. In particular, it 
is easy to check that the equation for pf^^z) j-gjj^ains the 
same as in the unshifted case, Eq. ([3]). 

We consider the exact calculation of the threshold \th 
for the existence of an active-to-absorbing phase transi- 
tion in BARW-DP. In order to check for the presence of 
a phase transition it is enough to study the behavior of 
A = r^-'^'^^ ij) = 0) as a function of the annihilation rate A, 
since the zeroes of A correspond to a divergence of the 
correlation length. Given that \th corresponds to the 
transition value of A when a — >■ 0'*', an analysis at lead- 
ing order in a allows for an exact calculation of Xth- An 
equation for T^^'^\p) at order Oia) can be represented 
in the diagrammatic form of Fig. SI which can be written 

T^'-^\p) = -a + Aal{p) I G{q)G{p - q) + 0{a^) (10) 

Jq 

where we have evaluated the propagator G{p) and the 
vertex r"'^-'(g,p — q,—p) at order zero in a, that is in 
PA, and consequently replaced this last function by 4:1 (p) 




FIG. 4. Closed equation for F'^'^) in BARW-DP at 0{a). 

By substituting the expression for l{p — 0), evaluating 
[T0| at p = 0, and putting A = 0, we find 

Xth ^ (2 1 G{q)G{-q)\ . (11) 



Xth/Da"^-'^ (this work) 3.96 6.45 8.65 10.7 

Xth/Da'^-'^ (MonteCarlo) [13| 3.99 6.48 8.6 10.8 

TABLE I. Values of the threshold coupling Xth for various 
dimensions d. 



To evaluate Xth, we need to take into account that the 
properties of a phase diagram are not universal, and de- 
pend on the specific form of the theory at short distances. 
This is as in equilibrium statistical mechanics, where crit- 
ical temperatures depend on the specific form of the lat- 
tice. Here we consider a particular microscopic form for 
the model corresponding to a hypercubic lattice with lat- 
tice spacing a. In this case the propagator reads 



G{q) = 



1 



#^Er=i(l-cos(ag,)) 



(12) 



In Table [H the value of the resulting threshold coupling 
Xth is given, proving in particular the existence of a phase 
transition in every dimension. Previous results from 
Monte-Carlo simulations 13|, |22| are in excellent agree- 
ment with these exact ones. This same general structure 
of the phase diagram has been also been show to exist in 
other models in the DP universality class [23| . 

It is convenient to point out that for d < 2 an IR di- 
vergence of the integral in (fTTj) takes place. This makes 
Xth = in those dimensions, in agreement with the re- 
sults of 2]. For this reason, for d < 2 it is not useful 
to expand the model at small a for a finite A in order 
to study the phase transition. Moreover, this also shows 
that in those dimensions the transition is dominated by 
IR effects, and correspondingly most of the dependence 
on the microscopic behaviour of the model is absent. 

Concluding remarks. In this work we take advantage 
of the structure of PA to find closed exact expressions for 
any of its response functions, which we use to perform 
an expansion in the branching rate a around this model. 
This gives us access to the small branching regime of 
BARW in two important universality classes. 

In the case of the system of reactions 2A — 5- 0, A — > 2A, 
which belongs to the DP universality class, we give an 
explicit proof of the existence of a phase transition in all 
space dimensions. We have moreover calculated exactly 
the non-universal threshold value of the annihilation rate 
for this phase transition to occur. This result is beyond 
the possibilities of usual perturbation theory. 

As for the parity conserving case, we find, surprisingly, 
that the appearance of the PC fixed point associated with 
an active-to-absorbing transition below a critical dimen- 
sion must occur at a nonzero value of the branching rate, 
which would be compatible with a scenario where there 
are not one but two new fixed points as d is lowered. 
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